Analytical solution for optimal squeezing of wave packet of a trapped quantum particle 
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Optimal control problem with a goal to squeeze wave packet of a trapped quantum particle 
is considered and solved analytically using adiabatic approximation. The analytical solution that 
drives the particle into a highly localized final state is presented for a case of an infinite well trapping 
potential. The presented solution may be applied to increase the resolution of atom lithography. 

OO I The recent interest in squeezed quantum states has stimulated research of optimal squeezing of quantum wave 

^^ . packets. Several schemes, including numerical solution of optimal control problem [l|, I4I, or iterative pump-dump 

~[^ ■ technique by switching the wavefunction between the ground and excited states with it pulses t3|, were introduced. 

Another method, based on classical parametric squeezing by a sudden increase of the potential depth was considered 

O ' for the case of 3D optical lattice gases [J]. Other, more efficient methods, need relatively short (5- like control pulses 

.^ • [5| . An experimental realization of this squeezing algorithm using impulse periodic potential |6| demonstrated the 

fundamental limitation for squeezing due to the nonzero duration of the kicking pulses. An interesting analytical 

^^ . solution was proposed for controls on time scales shorter than classical Kepler orbit period [TJ. However, most of the 

above mentioned methods were based on numerical solutions and have demonstrated limited success. 

1— I . One of the significant applications of squeezed atomic wave packets is atom lithography [8| . Atom lithography is 

1-5 ', relatively cheap and powerful method that does not require any material mask to fabricate semiconductor, metallic 

or magnetic nanostructures. The method is based on manipulation of atoms by laser fields on a scale under 100 nm, 

and it is considered as a promising way to overcome the resolution limits of usual optical lithography [5| . If we know 

J3 that at a moment T a cold atom will arrive at a surface, then the uncertainty of its position in the surface's plane 

^ may be roughly described in terms of the corresponding width of the atomic wave packet (we neglect diffusion on the 

O"^. surface and other effects that also have contribution to the uncertainty of the final atomic position) . 

Another possible application of squeezed molecular wave packets is mapping of the potential surfaces [Sj. Using 
pump-probe" measurement technique the excited state dynamics can be followed and the potential surface can be 
^ ' mapped. The ability to map potential surfaces using such techniques requires generation of very localized wave packets 

In this paper we are going to consider the problem of squeezing of a wave packet of a trapped quantum particle. 

X^ . Using our strategy one can obtain, in principle, an arbitrary narrow wave packet, that, in return, results, for example, 

in ultrahigh resolution of atomic lithography. However, it is clear that the squeezing of wave packets is limited by the 

Heisenberg principle. For example, a Gaussian- like wave packet with the width a has the kinetic energy Ekm oc a~'^. 

QQ Thus, one needs to provide the controlled system with enough amount of energy. The amplitude of the control field 

(^ : should be relatively big in order to achieve wave packet localization. As a result, perturbation theory may not be 

t^ ' applicable, and numerical solution of the time-dependent problem may be necessary. However, we are going to show 

that under curtain approximations the optimal control field and the corresponding wave packet evolution can be 

described analytically. 

3 . Let us consider a simplified ID picture, where a quantum particle of mass m is trapped in the potential Utrap{x) 

■ ■ ' and interacts with a time-dependent control potential Uc{x,t). Since in this work optimal control is achieved using 

linear polarized field and therefore can be performed in x and y directions independently, 2D generalization of the 

presented approach is straightforward. In the absence of Uc{x, t), the eigenstates and the eigenenergies of the trapped 

particle are 4'jix) and Ej (j = 1,2, ...) correspondingly. 

The dynamics of the particle's wave packet is described by the time-dependent Schrodinger equation 

ih—^{Ho + U,{x,t))^, (1) 
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where Hq = 2^-;^;^ + Utrap{x) is the unperturbed Hamiltonian. We assume that initially system is in the ground state 
^(a;,0) = il^i(x). Note, we neglect in this work any kind of decoherence effects (including spontaneous emission). If 
T is the duration of the control time interval and 7 is the characteristic decoherence rate in the system, we assume 



that 7T <^ 1. We consider the control potential in the form: 

Uc{x, t) = U^n{x) Y, y^t) cos(wijt), (2) 

where ljij = {Ej — Ei) / h is the transition frequency between levels 1 and j, and Vj{t) is the envelop of the component 
of the control field with the carrier frequency luij . Our guess of the time dependent structure of the optimal control 
potential Uc{x,t) is based on a natural assumption that the most efficient control is achieved within the resonance 
coupling ill| . The function Uin{x) characterizes the inhomogeneity of the control potential. As an example, one can 
consider form Uin{x) oc sm((3x)/(3 one and Uin{x) oc x represents the long wavelength limit /3 — * (linear potential, 
constant force). This potential can be realized in different ways. For example, one can use locally inhomogeneous 
electric field, so the force acting on a neutral atom is proportional to the gradient of the field intensity |dj}] . Another 
approach is to use ions instead of neutral atoms to control. In this case the force acting on an ion is simply proportional 
to its charge and the local amplitude of the applied field. 

The optimal control problem is usually formulated as the following: starting from the initial state ipi{x), one would 
like to steer the system optimally to a target state ^T(a;) = "^{x^T) at time t = T, such that the spatial dispersion 
a of the wave packet is minimal: a^ =< vE'tIx^JvE't > — < '^t\x\'^t >^^ min. Obviously, this condition on ^t is 
vague, because there is an infinite amount of wave packets with the same cr, and as a result, the numerical algorithm 
can be easily trapped in local minima. In other words, the condition a — > min allows very large search space, that 
may complicate iterative numerical solution of the optimal control problem (see for example, 1]). 

A simple way to handle this uncertainty is to determine a target wave packet with a specified shape, that can be 
parameterized by a. Thus, one can formulate optimal control problem that has a unique solution for arbitrary small 
a. The problem is then split into two parts, and the first part is to find an appropriate shape of the target wave 
packet. A reasonable choice can be made if we look for a strongly squeezed wave packet with a minimum kinetic 
energy for a given width. In this case we minimize the necessary work to be done to build up such a wave packet. It 
is easy to show, that the ground state of a harmonic potential satisfies this criterium. However, as we show below, in 
the limit of long wavelength control field the initial and the target wave functions should have different symmetry, i.e. 
if we start from the ground state, which is obviously a symmetric (even) function, one needs an antisymmetric (odd) 
target wave function. This is the reason to choose the target wave function in the form of the first excited state of a 
harmonic potential. 

We search for a control field that is optimal in the sense of spending minimum amount of energy that is still 
enough to reach the target state. Thus, we introduce a constraint that effectively limits the total field energy 
Si=^ /„ \Vj{t)\^dt — Etot- Note, that without this constraint optimal control problem has no finite solution, because 
as we mentioned above, one needs an infinite amount of energy Etot -^ +00 to squeeze a broad wave packet in order 
to achieve cr — > 0. 

Using expansion in a real eigenbasis of the trapping potential ^(x, t) ~ Tlij^i ^j(J') exp{—iEjt/h)ipj{x), Eq.([T]) can 
be written as 
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ihcij ^2^Vp{t)cos{uJipt)2_^ak{t)e * " dkj, (3) 

p=2 k=l 

where dkj is the coupling matrix element between states k and j: d^j — J_ iljk{x)Uin{x)'ipj{x)dx. 

The optimal control problem is to find the time dependent field envelopes Vj{t) that drive the system to the final 
state characterized_by the set of amplitudes {aj} at time T, aj — J_^ ipj{x)'^T{x)dx. The Lagrangian of this 
optimal control is p^j : 

L= fY.^\aAt)-a^?5{t-T) + \V^{t))dt, (4) 

where A is a Lagrange multiplier and 5{t) is the Dirac delta function. The problem of optimization of the Lagrangian 
Eq.Q subject to non-holonomic constraints given by Eq.Q is very complicated, and needs some simplifications. 
Let us assume the adiabatically slow changes of the field envelopes Vj{t) on the characteristic time scale maxjw^^}. 



Applying the Rotating Wave Approximation (RWA)[l3| to Eq.Q one gets: 

*^"J " 2 ^ Vp{t){A{nujip -Ek+ Ej) 

+A{-nujip-Ek + Ej))ak{t)dkj, (5) 



where the function A() returns unity if its argument equals to zero, and zero in all other cases. Thus, all the off- 
resonance quickly oscillating terms are set to zero. In Eq.® we account all resonance transitions, however, as we see 
later, we can neglect transitions between initially unoccupied states. 

Now we make a crucial approximation that significantly reduces the complexity of the problem, and it is based on 
the assumption that the target wave packet is narrow. If a is small, a/L <C 1 (L is the characteristic length scale of 
the trapping potential Utrap{x)), it is easy to show that the target amplitudes aj are also small, |aj| oc \/ cr jE ^ 1. 
It is natural that the overlap between an eigenfunction ^j and a narrow target wave packet is small. 

Since under the optimal control the amplitude of an excited state aj{t),j > 1 should grow monotonously from its 
zero initial value to its final value \aj\ ^ 1, one can simply neglect aj{t),j > 1 compare with ai(t). Thus, Eq.(I5]) can 
be simplified, since the |ai(i)| 3> |aj(t)| most of the time, except may be for a short time interval close to T. After 
these simplifications Eq.(l5]) becomes: 



1 ^ 
*^"i " 9 X! ^'^ (^^^'^ it)dki , (6) 



fc=2 



ihcik = -Vkit)aiit)dki. 

We have performed comparison between numerical solutions of Eq.((5]) and Eq.([6|), and for small parameter values 
(7 < lO^^L found a good agreement (< 1%). 



Using vector representation of Eq.([5|): a — Z{t)a and the Magnus series expansion 1J|, we can write the exact 
solution as 

a(t) = e L J 0,(0), (7) 

where vector a(0) determines the initial wave packet "^{x, 0), and [ , ] denotes the commutator. 

We neglect all the terms in the series Eq.© except the first one. It is a good approximation, if the first term in 
the series is much larger than the second one, and all further terms. Using the explicit form of the operator Z{t) and 
the Mean Value Theorem, we obtain |^(0| ^ ^^| — dt\ ' '^-'^^re t,ti 6 (0, T). This is the adiabatical condition, 
under which it is possible to integrate Eq.® analytically. Introducing new variables Ok{t) — 1/2 /„ Vk{t')dt', and 

-^(0 = ( J2k^ ^ki^)'^ki ) ' O'^s can write an approximate solution given by Eq.([7]) as: 

ai(i) ^ cos (R{t)/h\ 
aj{t) = d jiO j (t) sin (R{t)/fi\/R{t) J = 2,3,... (8) 

We can use this solution to substitute into the Lagrangian Eq.(j4]), that allows us to derive explicit Euler-Lagrange 
equations with respect to unknown Ok(t). Note, that the approximations made in Eqs. (|5l7p do not affect the normal- 
ization condition: X]fe=?i l^*:(OP — 1j unlike if one uses the simple first order perturbation theory. 
The Euler-Lagrange equations which determine optimal 0fc(i), are: 

^^Ok - ^ E l«^({^^-}) - «Jl''^(^ -T) = 0, (9) 

j 

where we use the explicit notation aj({6k}) to stress the dependence of amplitudes on the set of unknown functions 
{Ok}- 

Since in practice the target state has a finite (although may be very small) width a, it is reasonable to keep a 
limited number of energy levels N of the controlled system in the consideration. For a smaller a one needs to increase 
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FIG. 1: (Color online). Time evolution of the controlled wave packet density |^(x,t)p in an infinite rectangular well. Note monotonic 
increase of the wave packet height and decrease of its width. 

N. The appropriate number N can be determined from the condition J2i=^+i kf P < Cconst «C 1. 

Note, that under the discussed approximations optimal control problem Eq. ([9]) in new variables becomes particular 
simple and formally splits into iV — 1 independent control problems for the whole control interval (since the delta 
function takes zero value), except at the time t ~ T . Each of the equations is equivalent to find an optimal envelope 
Vk(t) for one resonant transition. 

Integration of Eq.® gives Ok ^ Ak + Bkt. This result is obtained under the condition that A 7^ 0. Thus, linear 
dependence of 9k on time is a consequence of the constraint on the energy, introduced in Eq. ([4]) . Ak can be readily 
determined from the initial condition 6k{0) — 0, that gives Ak — 0. The coefficients Bk are determined from the 
condition that Eq. ^ should also be satisfied at t = T. This is equivalent that the second term in Eq. ([9|) should 
turn exactly to zero at t = T. This gives us a system of algebraic equations for Bk'- 

a,{{9k})\t=T = d.iB^T sin (^R{T)/h^/R{T) = aj,j = 2,3,... (10) 

/ \l/2 

with R{T) — Ti Y^k^ ^k^ki ) ■ Eq. (fTO|) actually plays a role of the second boundary condition for each dk(t) at 
t — T. It is easy to verify that the solution (see solution of a similar problem in [12]) 

Okit) = Bk = IVkit) = ^,k =. 2, 3, ... (11) 

satisfies Eq. (|10p . Here we have assumed that the symmetry of the target wave packet is chosen such, that if dki — 
then a[ = 0, and the correspondent Vfc(t) is also zero. 

The result Eq. fTTI) self-consistently justifies the adiabatical approximation made for the Magnus expansion Eq.©, 
since for the constant field envelopes all terms except the first one in Eq.© are equal to zero. Substituting Eq. (fTT|) 
into Eq. Q one easily gets the corresponding dynamics for the occupation numbers ak (t) . Using the obtained solution 
one also can estimate that the approximation made for Eq.® is a good one for times i G [0,T — 2T \Ja jLjii). 

The cost of all the simplifying assumptions (including RWA), which lead to Eq. ([9]), is that Eq. (I9|) is a system of 
the second order differential equations for unknown 0^, and one can not impose more than two boundary conditions 
for each equation. The constraint on the energy can be seen as a third condition, but for the whole system. By 

substituting into the energy constraint X]fc°=2 /n l^fcl^*^^ — Etot, one obtains L J2T=2 l^fcl^*^^ ~ '^'k'=2 l^fel^^ — Etot- 
Since Bk are already determined from the condition at t = T, the energy constraint Etot cannot be arbitrary. However, 
one can assume that duration of the control T is unknown at the beginning, so one can determine the control time T 
through a given total energy of the control field Etot ■ 

Now let us consider an example of a trapped quantum particle in an infinite potential well of width L, interacting 
with a control potential Uc{x, t) in the long wavelength limit /3 — > 0. As we mentioned above, it is easy to show that 
the matrix element dji —< V'fc|C^m|V'i >= ^4i(cos(7rj)j + j)/{Tr^{l — 2p + j^))^ so dji — 0, ii j is odd, and dji ^ 
if j is even. We assume the system initially is in the ground state, thus it is impossible to drive the system into a 
squeezed target state with the same symmetry in the long wavelength limit. In order to utilize the only allowed ground 
state-odd states transitions, the target state ^^rix) is chosen to be an antisymmetric wave packet parameterized by 



a and located at position x^: 

■^t{x) = B{x - xo) exp(-(a; - x^f/a'^), (12) 

where the constant B is determined from the normahzation condition over the interval [0, L]. In the limit cr ^ 0, 
'^xix) becomes equivalent to the first derivative of the Dirac delta function i5'(a; — xg), and it has the minimum kinetic 
energy among all other wave packets with the same symmetry and given width a. The amplitudes of the final state 
in the basis of the infinite well potential a^ can be calculated analytically in terms of the Erf(x) function. 

The ultimate goal of a useful squeezing algorithm is to create a wave packet with a characteristic size much less than 
the characteristic length scale of the inhomogeneous control potential. A possible realization of the proposed squeezing 
method in the long wave length limit is to control ions of charge +q in the oscillating linear potentials. Let us assume 
for simplicity g = |e|. In a typical experiment let us assume the intensity of the highest mode of the controlling 
electromagnetic field is of the order of 100 mW/cm^, that corresponds to the electric field of ii^ = 8.0 x 10~^V/m. We 
assume the initial width of the wave packet to be of the order of the infinite well trapping potential L = 1 x 10~^m, 
and the final width of the wave packet a k, L/50 = 20 nm. Using Eq. pTJ) . we estimate duration of the squeezing 
procedure T = Trha'^Q/{Edi 30) « 3 x 10""' s. 

In Fig. 1 we show time evolution of the wave packet density |^(x, i)p, trapped in the infinite well potential x e (0, L) 
under the control of the optimal field determined by Eg. pTjl . We choose A^ = 30 of controlled levels, xq = L/2, and 
final width a = L/50. Note, that the height of the controlled wave packet is increasing monotonically, while the width 
is monotonically decreasing to its theoretical target limit. This is a direct consequence that the control is optimal, 
opposite to the previous studies (see, for example, Fig. 3 in [l|, the non-monotonous decrease of the width of the 
wave packet is an indication that the numerical solution is a local extremum). 

If it is more preferable to use neutral atoms instead of ions, then the necessary control potential can be created 
using dipole forces which are proportional to the amplitude gradient of the laser field [10|. Note, that unlike in the 
case of ions, the dipole force arises for sufficiently inhomogeneous laser fields. 

Another possible way to realize the proposed squeezing method is to consider control of molecular wave packets, or 
electron wave packets in artificial molecules: quantum dots. In this case one starts from the ground state and finishes 
at the linear combination of the vibrational (or eigenstates of the confining potential in the case of quantum dots) 
states, that builds up a narrow wave packet at a given time T. The latter can be probed by means of the ultrafast 
spectroscopy. 

In this work we have derived analytically the control potential to localize an atomic wave packet in real space. This 
method is based on our knowledge of the eigenenergies and eigenfunctions of the controlled system in the trapping 
potential, and it is rigorous under the approximations made. Note, we can achieve the wave packet squeezing using 
control potentials with much larger wavelength (/3 — > 0) compare to the characteristic scale of the squeezed wave 
packet a, that may be useful in atom lithography. 

One essential assumption we made is a strong unharmonicity of the trapping potential. This assumption assures 
us that resonance frequencies between the initial (ground) state and multiple final states are all well separated. The 
trapping potential can be different from the infinite well potential considered in this work, but it will be more difficult 
to obtain solution in a simple analytical form. The presented method breaks down for relatively short control intervals 
T or relatively strong control fields, when the KWA is not applicable. Another limitation is that the control interval 
cannot be very large, since the approximation that decoherence effects are negligible T-f -C 1 should hold. 

The resulting optimal control field has a complicated spectrum and shape, and it may be hard to obtain it using the 
standard numerical solution techniques for optimal control problems [l'| . That makes analytical form of the presented 
solutions even more attractive. 

This work was carried out under the auspices of the National Nuclear Security Administration of the U.S. Depart- 
ment of Energy at Los Alamos National Laboratory under Contract No. DE-AC52-06NA25396. 
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